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Characteristics of the muon component in EAS are analyzed together with their fluctuations. The aim of this 
analysis — a comparison of experimental data with computational results obtained within frameworks of various 
hadron interaction models for protons and iron nuclei and an estimation of cosmic ray mass composition in the 
ultra-high energy region. 



1. Introduction 

The Yakutsk complex array for many years mea- 
sures three main observables of extensive air show- 
ers (EAS): total charged component, muons with 
£thr > 1 GeV and Cherenkov light [lj. Using these 
data we estimated the EAS energy with model inde- 
pendent quasi-calorimetric method [1| and determined 
the depth of maximum shower development (by the 
measured Cherenkov light lateral distribution, using 
the parameter p — lg Q(200)/Q(550) and by the shape 
of the Cherenkov light pulse, t 1 / 2 ) 0, E|- The relative 
muon content at different core distances was measured 
d, Q and the cosmic ray (CR) mass composition was 
estimated by various EAS characteristics 

In this paper we analyze the muon component of 
EAS: mean characteristics, muon content and its fluc- 
tuations at fixed energy. The analysis is conducted 
within the framework of QGS Jet II and EPOS [TO] 
hadron interaction models involving computations for 
primary particles of different masses using CORSIKA- 
6.900 code [12. 



2. Muon lateral distribution function 

In Figure [T] examples of mean lateral distributions 
for muons at different energies are displayed. The 
muon lateral distribution function (LDF) is signif- 
icantly lower then that of the charged component 
and can be effectively measured in individual events 
at Eq > 10 17 eV within the core distance range 
100—800 m. Thus, as a classification parameter in this 
energy region, a parameter p M (600) could be used — 
the density of muon flux at 600 m from shower core. 

Solid and dotted lines on the figure denote computa- 
tional results obtained with QGSJct(UrQMD) models 
for proton (solid) and iron (dotted). It is seen from 
Figure Q] that the muon LDF from protons is steeper 
than that from iron nuclei and this difference is es- 
pecially pronounced at large core distances. Quali- 
tative comparison of computational results with the 
experiment reveals a better agreement with a heavier 
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Figure 1: Lateral distribution of charged particles for 
fixed energy values 10 17 , 10 28 , 10 19 eV and zenith angle 
9 = 15°. Symbols refer to the Yakutsk experiment. Solid 
lines denote computational results with 
QGSJet+FLUKA [HI for protons, dotted line — for 
iron, red line — EPOS+UrQMD P3 for proton, green — 
carbon, blue — iron. 

component of primary CR at Eq < 10 18 eV and with 
lighter at Eq ~ 10 19 eV. This feature could be stressed 
out if one puts parameter r 2 ■ p(r) on the y-axis of a 
plot instead of simple p{r). 

3. Muon portion and its dependence on 
angle, energy and the depth of maximum 
EAS development 

We considered the dependence of p^j p c h on the 
length of shower development after the maximum — 
A A = xo/cosf? — i mM . In highly inclined showers 
the muon content increases proportionally to xq / cos 9 
value, where xq = 1020 g/cm 2 for Yakutsk. 

It is a known fact that the depth of maximum EAS 
development differs significantly, depending on the 
kind of primary particle and, hence, this feature could 
be used in the analysis of the CR mass composition: 
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for instance, by fixing the AA parameter and study- 
ing the fluctuations of p^j p c n value. This technique 
is rather similar to one proposed by Atrashkevich et 
al [H|. 

Shower parameters calculated with CORSIKA code 
were modified by applying distortions according to ex- 
perimental errors. Parameters measured in experi- 
ment (e.g. cos0,2: roaX) p c h(r), p M (r)) for every shower 
were rolled with the normal distribution with a pa- 
rameter according to the experiment: 

cr(0) =3 -sec (9; 
o-(xmax) = 40g/cm 2 ; 

*(pr) = ( 0-025 + — -) 

V V s dot ' Pr ' COS 9 J 

where Sdet is the area of the detector. 

Figure [5] shows the dependence of p^/pch on the 
length of cascade development after the shower maxi- 
mum compared with computational results. A strong 
correlation is observed between the muon content and 
the length of track in the atmosphere after the shower 
maximum. It is also seen that experimental data are 
in good agreement with simulation results. 



4. Mean characteristics 

Figure shows the energy dependence of p p (600) 
obtained in the Yakutsk and MIA EAS experiments 
16]. A good agreement is observed between the two. 

Computational results obtained with EPOS and 
QGSJctOl from the work by Abu-Zayyad et al 
are denoted with lines, dotted line represent our sim- 
ulations with QGSJct II for protons and iron nuclei. 
A comparison of our computations with the results 
obtained by Abu-Zayyad et al [16[ reveals that vir- 
tually there is no difference between QGSJetOl and 
QGSJet II. A significant discrepancy is observed be- 
tween EPOS and QGSJet II and it, as we believe, is 
associated with different amounts of muons generated 
by models. For example, the /9 M (600) value calculated 
with EPOS for proton coincides with ^(600) obtained 
with QGSJet II for iron. Thus, a comparison of exper- 
imental data with model calculations result in contro- 
versial conclusions on CR mass composition. Accord- 
ing to EPOS, at energies up to 2 x 10 17 eV CRs consist 
of iron nuclei and above that energy, up to 10 19 eV - 
of protons. In the energy interval 10 17 — 3 x 10 18 eV 
QGSJet II computations agree with the experiment 
quite well if primary particles are iron nuclei and 
above 3 x 10 18 eV the mass composition might be 
mixed with portion of protons and helium nuclei not 
less than 50 — 60%. More precise estimation of CR 
mass composition could be derived after improvement 



of theoretical models and selecting one, that describes 
experimental EAS data better then others. 



5. Fluctuations of p^/p C h relation on the 
ground level at the energy ~ 10 18 eV 

Showers initiated by different nuclei have differing 
altitudes of the maximum which in turn means that 
different numbers of muons are generated in these 
showers. It also means that they cover different paths 
in the atmosphere. By analyzing the tracks of muons 
that they pass in the atmosphere after the maxi- 
mum of shower development we can try to estimate 
the composition of cosmic rays. With this aim in 
view, by choosing the mean zenith angle 36° (which 
corresponds to the track length after the maximum 
AA = 500 g/cm 2 ), let us normalize the values of muon 
content to this level and consider their fluctuations. 

Results are presented in Figure [4] Also shown are 
the computational results obtained with QGSJet II 
and EPOS models for various nuclei. Measured val- 
ues of fluctuations are presented in Table [I] Obtained 
results have shown that within this method fluctua- 
tions of p A1 (600)/p c ; l (600) parameters do not allow to 
estimate the CR mass composition. However, mean 
values from different nuclei differ. Besides, QGSJet II 
hints at a heavier composition than that of EPOS: ac- 
cording to first one, the composition of selected show- 
ers shifts towards heavier nuclei; according to sec- 
ond one, showers correspond to nuclei of intermediate 
group. On the whole, both models argue for a mixed 
composition. 

However, if one takes into account gamma-photons 
generated in ground covering muon detectors, the 
mean value of p M (600)/p c ^(600) relation decreases and 
the composition shifts towards lighter nuclei (protons- 
helium-carbon) [l8j . 



6. Conclusions 

Within the framework of QGSJet II and EPOS 
hadron interaction models using the CORSIKA code 
the values of muon portion Pp,/p c h at core distance 
r = 600 m were obtained. A relation between the 
muon portion and a distance to the depth of shower 
maximum AA was also obtained. A comparison of the 
dependency with experiment has shown that taking 
account of the experimental errors in the simulation 
data, a good agreement is observed between simula- 
tion and experiment. 

A comparison of the muon portion distribution with 
computational results points towards a mixed cosmic 
ray composition near E$ > 10 18 eV. Large fluctua- 
tions of the muon portion prevent revealing of a single 
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Figure 2: A dependence of muon portion with e t hr > 1 GeV on the length of track in the atmosphere for individual 
showers with 8 — — 50° and energy 10 18 eV. On the left — results obtained with QGSJet II model, on the right — 
with EPOS model. 



Table I Fluctuations of p M / p c h relation, normalized to the length of track 500 g/cm 2 

QGSJet II, FLUKA EPOS, UrQMD 

Yakutsk Yakutsk" p C Fe p C Fe 

(jWPch) 0-3145 0.2768 0.2687 0.3025 0.3193 0.2893 0.3170 0.3381 

a 0.0747 0.0657 0.0517 0.0541 0.0536 0.0563 0.0511 0.0539 

"With respect to contribution from gammas generated in the 
shielding of detector (Dedenko, 2010) 
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Figure 3: Muon density p M (600) measured in Yakutsk 
experiment as a function of primary energy compared to 
model calculations. 



determined group of nuclei. A more detailed analy- 
sis is required, involving possible systematics of muon 
density measurement in the Yakutsk experiment. 
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Figure 4: A distribution of the p M /p c h relation, normalized to the track length 500 g/cm 2 . On the left — according to 
models QGSJet II(FLUKA), on the right — according to EPOS(UrQMD). 
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